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Abstract. The Hipparcos Transit Data are a collection of 
partially reduced, fully calibrated observations of (mostly) 
double and multiple stars obtained with the ESA Hippar- 
cos astrometry satellite. The data are publicly available, 
as part of the CD-ROM set distributed with the Hip- 
parcos and Tycho Catalogues (ESA SP-1200, |1997|) , for 



about a third of the Hipparcos Catalogue entries includ- 
ing all confirmed or suspected non-single stars. The Tran- 
sit Data consist of signal modulation parameters derived 
from the individual transits of the targets across the Hip- 
parcos focal grid. The Transit Data permit re-reduction of 
the satellite data for individual objects, using arbitrarily 
complex object models in which time-variable photomet- 
ric as well as geometric characteristics may be taken into 
account. We describe the structure and contents of the 
Transit Data files and give examples of how the data can 
be used. Some of the applications use standard astronomi- 
cal software: Difmap or AIPS for aperture synthesis imag- 
ing, and GaussFit for detailed model fitting. Fortran code 
converting the data into formats suitable for these applica- 
tion programs has been made public in order to encourage 
and facilitate the use of Hipparcos Transit Data. 

Key words: Methods: data analysis - Space vehicles - 
Techniques: image processing - Catalogs - Astrometry - 
Binaries: general 



1. Introduction 

In the course of its three-year lifetime the European Space 
Agency's Hipparcos satellite performed some 13 million 
scans across the ~ 118 000 stellar objects on its pre- 
defined observing list. For the vast majority of targets 
the subsequent data reductions succeeded in determining, 
accurately and without ambiguity, a small set of astromet- 
ric and photometric parameters which contain practically 
all the useful information on the objects. In the simplest 
case of a non-variable single star these data would be the 



position of the star at a certain epoch, the trigonomet- 
ric parallax, the two proper motion components, a mean 
magnitude, and various statistics such as the mean errors 
of all the parameters. For thousands of variable, double 
and multiple stars, the necessary additional parameters 
were determined without too much problem. All of these 
results are readily available in the Hipparcos and Tycho 
Catalogues (ESA [l99^ ). For most objects and for most ap- 
plications of the results, there is no need to probe deeper 
into the intricacies of the Hipparcos data acquisition and 
processing. 

It was inevitable, however, that the observations of 
some (very small) fraction of the objects could not eas- 
ily be interpreted in terms of standard models, and that 
inadequate, or even erroneous, object models were ap- 
plied in some cases. Potentially very valuable information 
could have been lost for these objects if only summary re- 
sults, based on inadequate models, were published. This 
situation has been avoided by including a good selection 
of intermediate results in the published catalogue. Being 
only partially reduced these results are less dependent on 
specific model assumptions, but still relatively simple to 
use thanks to a careful calibration into the photometric 
and astrometric systems of the final catalogue. The pub- 
lished version of the Hipparcos and Tycho Catalogues in- 
cludes six CD-ROMs, five of which contain mainly such 
calibrated, partially reduced data. The data fall in three 
distinct categories: 

— purely photometric information: the Epoch Photome- 
try, consisting of calibrated magnitudes for each tran- 
sit of the targets across the field, allowing a posteriori 
analysis e.g. of stellar variability (van Leeuwen et al. 



1997) 
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* Based on observations made with the ESA Hipparcos as- 
trometry satellite 



purely geometric information: the Intermediate Astro- 
metric Data, which are the one-dimensional coordi- 
nates ('abscissae') of the targets on reference great cir- 
cles, each coordinate derived from several consecutive 
transits and referring to a (Hipparcos-specific) centroid 
of the target (van Leeuwen & Evans 1998| ) ; 
combined photometric and geometric information: the 
Transit Data (TD), which are the Fourier coefficients 
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(equivalent to amplitude and phase) of the light mod- 
ulation of the individual transits for each of the targets 
across the 0?9 x 0?9 main grid of the Hipparcos instru- 
ment. 

These mission products are fully described in Sects. 2.5, 
2.8 and 2.9 of Vol. 1 of the Hipparcos and Tycho Cat- 
alogues. Examples using the Intermediate Astrometric 
Data have been given by van Leeuwen & Evans (1998). 

While the TD thus provide information that is more 
detailed and on a more fundamental level (in the sense 
of being closer to the satellite 'raw' data) than either the 
Epoch Photometry or the Intermediate Astrometric Data, 
two important restrictions should be noted: (1) TD are 
only available for about a third of the targets, or some 
38 000 objects; (2) TD are based solely on the results from 
one of the data reduction consortia, the Nort hern Data 
Analysis Consortium (NDAC; Lindegren et al. 1992 ). By 
contrast, the Epoch Photometry and the Intermediate As- 
trometric Data are available for virtually all the 118 204 
catalogue entries, and results from both NDAC and the 
Fundamental Astrometry by Space Techniques Consor- 
tium (FAST; Kovalevsky et al. 1992) were normally used. 

Detailed satellite operations and data reductions are 
described elsewhere (e.g., Perryman & Hassan 1989] ; Per- 
ryman et al 
al 



1989 



Perryman et al 
19921; Lindegren et al. |1992|; Vols 



1992; Kovalevsky et 



2-4 of ESA 1997 



van Leeuwen 1997). It will however prove useful to re- 



view how the TD were measured in order to understand 
their precise meaning (and limitations), aiding in the full 
exploitation of the information. Thus, Sect. 2 describes 
the convolution of the stellar diffraction image with the 
modulation grid, and the resulting detector signal mod- 
eled by the TD. Section 3 is a detailed description of the 
contents and format of the published TD files. Two ap- 
plications of the TD using publicly available software are 
demonstrated, viz. aperture synthesis imaging (Sect. 4) 
and model fitting (Sect. 5). 

2. What are the Transit Data? 

2.1. Availability of the data 

The TD are contained on a CD-ROM (Disk 6) in Vol. 17 
of the Hipparcos and Tycho Catalogues. A formal descrip- 
tion of the TD, including detailed format specifications of 
the CD-ROM files, is found in Vol. 1, Sect. 2.9, of the 
Hipparcos and Tycho Catalogues. The TD come from an 
intermediate step of the data reductions performed by 
NDAC. In previous publications, the equivalent interme- 
diate NDAC data were referred to as 'Case History Files' 
(Soderhjelm et al. 



1992). The data represent the scans 



(transits) of a selected number of targets from the Hip- 
parcos Catalogue (HIP), comprising 37 368 systems with 
38 535 different HIP entries. (Some systems have two or 
three HIP entries; cf. Fig. 0.) All Hipparcos stars that were 
classified as double, multiple or suspected non-single are 



included in the TD. Also, stars having problematic solu- 
tions in the Hipparcos Catalogue are included in the data 
set. In order to provide reference objects for calibrations or 
comparisons, several thousand (~ 5 000) 'bona fide' single 
stars were also included in the TD. 



2.2. The modulated detector signal 

The physical layout of the Hipparcos optical instrument 
was a telescope with two viewing directions in a plane 
perpendicular to the satellite's spin axis. The two differ- 
ent fields of view were combined onto a single focal surface 
by a special combining mirror. The light from each pro- 
gram star within either field of view was focused onto the 
modulation grid, which was located in the focal surface. 

The modulation grid consisted of a series of opaque and 
transparent bands. As the satellite spun, the diffraction 
image of each star traversed the grid perpendicular to the 
bands, resulting in a periodic (~ 7 ms period) modulation 
of the light intensity behind the grid. The varying intensity 
was measured by the Image Dissector Tube (IDT) , a pho- 
tomultiplier with an electronically steerable sensitive spot 
('Instantaneous Field of View', IFOV) of about 30 arcsec 
diameter. During each 'interlacing period' of ~ 130 ms the 
IFOV would cycle through the programme stars located 
within the 0?9 x 0?9 field of view. There were on average 4 
to 5 programme stars within the field of view at any given 
time. 

The telescope entrance pupil (for each of the two view- 
ing directions) was semi-circular, with diameter 0.29 m 
and with some central obscuration. The Airy radius of 
the (ideal) diffraction image was thus around 0.5 arcsec 
for an effective wavelength of 550 nm. The modulation 
grid had a basic period of 1.2074 arcsec, i.e. the separa- 
tion between the centres of adjacent transparent bands 
(slits). The slit width was about 0.46 arcsec, well matched 
to the Airy disk size and representing a compromise be- 
tween a sharp intensity maximum (requiring narrow slits) 
and high photon throughput (requiring wide slits). 

Since the grid was periodic, the detector signal, be- 
ing the convolution of the diffraction image with the grid 
transmittance, must also be a periodic function of the im- 
age centroid coordinate on the grid. Disregarding noise 
and variations in detector sensitivity, etc., the signal there- 
fore consisted of a constant (DC) component plus modu- 
lated components having spatial frequencies that are inte- 
ger multiples of the fundamental grid frequency. However, 
from the convolution theorem it follows that the signal 
cannot contain higher spatial frequencies than were al- 
ready in the diffraction image. The maximum frequency in 
the diffraction image is given by the maximum separation 
of any two points in the pupil and the minimum detected 
wavelength (~ 350 nm). For a pupil diameter of 0.29 m 
this gives a spatial period of 0.25 arcsec. Thus the theo- 
retically highest frequency in the signal is four cycles per 
grid period (the 'fourth harmonic'). The fourth and third 
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harmonics, with cycles of 0.3-0.4 arcsec, are however very 
strongly damped by the slit width (0.46 arcsec), which 
causes an averaging over more than one cycle for these 
components. As a result, the detector signal in practice 
contains only the first and second harmonics, in addition 
to the DC (mean intensity) component. Given the spa- 
tial frequency, the detector signal is therefore completely 
parametrized by five numbers, i.e. the DC term and the 
coefficients of the first two harmonics in the Fourier series 
representing the periodic signal. Furthermore, since Pois- 
son (photon) noise is by far the dominating noise source 
for most Hipparcos observations, it can be shown that a 
proper estimate of these five Fourier coefficients consti- 
tutes a sufficient statistic for the further estimation of the 
photometric and geometric characteristics of the target, 
independent of its complexity. The TD contain precisely 
these Fourier coefficients along with the spatial frequencies 
and other ancillary data. From the viewpoint of statistical 
estimation, practically no information was therefore lost 
by compressing the raw photon counts (on average some 
4500 bytes per transit and target) into the five Fourier 
coefficients (included in the TD file). 

The detector signal thus measured certain components 
of the diffraction image moving over the modulation grid. 
The diffraction image, of course, had aberrations due to 
the imperfections of the telescope and chromatic effects 
caused by the wavelength-dependent diffraction. Also, the 
sensitivity of the IDT varied across the field and over the 
time. All of these factors affected the detector signal. The 
user of the TD need not worry about this, because the 
TD have been 'rectified', which means all instrumental 
and colour effects have been removed as far as possible by 
means of the various calibrations produced in the data re- 
ductions. Therefore, the TD should represent the response 
of an idealized, constant, and well-defined instrument to 
the object. 

2.3. The reference point 

One of the key points of the TD is that the modulation 
phase in each scan is expressed with respect to a well- 
defined reference point on the sky. The astrometric pa- 
rameters for the reference point (different for each target) 
are also specified in the TD. Usually the reference point 
corresponds to the data gi ven in the Hipparcos Input Cat- 
alogue (HIC; Turon et al. 1992 ). However, the phase cal- 
ibration of the TD was made as if the coordinates of the 
reference point were expressed in the Hipparcos reference 
frame (nominally coinciding with the International Celes- 
tial Reference System, ICRS; Feissel & Mignard |199SP . 
The phase offset of the TD from the reference point thus 
gives the differential correction to the position from the 
reference point in the ICRS system. 

The reference point in general has its own values 
for proper motion and parallax, in addition to position. 
These values are contained in the TD 'header record', see 



1.2074 
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Fig. 1. This figure illustrates the definition of reference 
phase p (as in Eq. [j]) by means of the reference point R on 
the sky (with specified astrometric parameters), p = at 
the instant when a slit is exactly centred on the reference 
point. In the lower diagram the curves A and B show the 
detector signal, as function of reference phase, produced 
by the same point source S when scanned in different di- 
rections. The intensity maximum is displaced from p = 
depending on the separation of S and R as projected on 
the grid, modulo the grid period 1.2074 arcsec. 



Sect. 3.2.1. Of course, the real star also has a position, 



proper motion, and parallax. These two points move inde- 
pendently on the sky, according to their respective proper 
motion and parallax. In a particular transit the phase is 
determined by the scan direction and the distance between 
the two stars perpendicular to the slits. Figure [l] illustrates 
how the phase varies in two different scans across the star 
(S) and reference point (R). In scan A the star is in phase 
with the reference point, producing a light maximum for 
p = (and more generally for p — n2ir, where n is an 
integer) . Scan B has the star out of phase with the refer- 
ence point, resulting in a light maximum about a third of 
a period later (p ~ 2 rad). 

The raw detector signal consists of a sequence of pho- 
ton counts, TYjt, k — 1,2, ... , obtained in successive sam- 
ples of 1/1200 s integration time. The counts represent 
an underlying deterministic intensity modulation which 
(after background subtraction and various calibrations) is 
modeled as 

h = h + b 2 cosp k + b 3 sinp fe + 6 4 cos 2p k + 6 5 sin 2p k (1) 

where 61-65 are the Fourier coefficients given in the TD. 
Jfe is the expected stellar count rate in sample k (expressed 
in photons per sample of 1/1200 s) and p k is the reference 
phase of the sample. The reference phase is defined by the 
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following two conditions (cf. Fig. [l]): (1) that pk = if a 
slit is exactly centred on the reference point at the mid- 
time of sample k; and (2) that pk is increasing from to 
2tt over the time it takes the grid to move one grid period 
(1.2074 arcsec) over the reference point. 

Equation ([!]) is the most general representation of the 
signal produced by any object. In order to interpret this 
signal in terms of an object model it is necessary to know 
what the signal would be for a point source of unit in- 
tensity; in other words, we need the equivalent of the 
point-spread function i n im age processing. The rectifica- 
tion mentioned in Sect. 2J2 meant that the expected signal 
from a point source of magnitude Hp — (in the Hippar- 
cos photometric system) located at the reference point is: 



I k = K[l + Mi cosp fc + M 2 cos 2p k ] 



(2) 



Here K = 6200 counts per sample, M x = 0.7100 and 
M2 = 0.2485 are constants chosen not far from the actual 
mean calibration values at mid-mission (see Figs. 14.4, 



14.3 and 5.6 in Vol. 3 of ESA 1997) 



2-4- Spatial frequencies 

As illustrated by Fig. [l] the resulting signal will depend 
on the direction of scanning, even in the simple case of 
a point source with a fixed position relative to the refer- 
ence point. The direction of scanning across the reference 
point can be specified by a position angle 9 defined in 
the usual way, i.e. with 9 = for a scan where the grid 
is moving in the direction of increasing S ('North'), and 
9 = 90° for a scan in the direction of increasing a ('East'). 
The spatial frequency of the grid, / (expressed in radi- 
ans of modulation per radian on the sky), is nominally 
(1296000 arcsec)/ (1.2074 arcsec) ~ 1073 400. In reality it 
varies slightly, mainly because of differential stellar aberra- 
tion, which changes the apparent scale by up to ±0.01 per 
cent depending on the barycentric velocity vector of the 
satellite. For the production of the TD data, the appar- 
ent scale and orientation of the grid, as projection onto 
the sky in the vicinity of the reference point, were strictly 
calculated for each transit, taking into account aberration 
as well as the satellite attitude and the calibrated field-to- 
grid geometric transformation. The result is expressed as 
the rectangular components of the spatial-frequency vec- 
tor in the tangent plane of the sky: 



fx = -/sin( 



fy = ~/cOS< 



(3) 



The minus signs mean that the spatial-frequency vector 
(/) points in the opposite direction to the motion of the 
grid across the reference point. This is purely a matter of 
convention, which was adopted for historic reasons. 

2.5. The phase of an arbitrary point source 

Consider now a point source of magnitude Hp which, in 
a particular scan, is displaced by (x, y) from the reference 



point, where x is measured positive towards increasing a 
and y towards increasing S. Given the spatial frequency 
components of the scan, f x and f y , the expected signal is: 

h = lQ-° AHp K[l + Mi cos(p k + </>) + M 2 cos 2( Pk + <f>)} 

(4) 

where 

<t> = f x x + f y y (5) 

is the phase shift on the grid caused by the positional 
offset. [That is added to pk in Eq. (Q), rather than sub- 
tracted, is consistent with the sign convention adopted in 
the definition of f x and f y .] 

Equations (B)-(]5j) are based on a linearization of the 
transformation between spherical coordinates and their 
projection on the tangent plane of the sky through the 
reference point. Within angular radius p of the reference 
point, the neglected non- linear terms are generally of or- 
der p 3 , or < 0.1 mas for p < 160 arcsec. Since the TD 
for a given target are normally confined to an area set by 
the size of the IDT sensitive spot (~ 30 arcsec), or a few 
times this area (for multiple pointings), the linearization 
is always adequate. 

The positional offset (x, y) is in general caused by a 
time-dependent combination of the differences between 
the astrometric parameters of the star and of the refer- 
ence point. The differences in a, 8, p a * and ps are easily 
converted to an offset which is a linear function of time.[] 
The parallax difference produces a shift which is more 
complicated to calculate, as it depends on the position 
(s) of the satellite relative to the solar system barycentre 
at the time of observation. To facilitate this calculation, 
a third spatial frequency component, f p , is supplied with 
the TD. This is basically just the scalar product —s'f, 
where the satellite position s is in astronomical units and 
/ is the previously defined spatial-frequency vector of the 
grid. With this definition it is found that a parallax dif- 
ference of Att relative to the reference point causes the 
additional phase shift A(f> = f p Air in Eq. (||). The com- 
plete expression for the phase of a star with astrometric 
parameters (a, 8, 7T, p a *, p,g) is therefore 

= f x (Aa* + tAp a *) + f y {AS + tAns) + f P Air (6) 

where t is the time of the transit, expressed in (Julian) 
years from J1991.25, and Aa* = (a — ao)cos<5o, AS = 
8-6 , An = 7r-7ro, Ap a * = p a *-p a *o, Afis = ^6~^S0 are 
the differences with respect to the astrometric parameters 
of the reference point, (cto, 8q, 7r , p a *o, M<5o)- 

The periodicity of the grid means that, in a given scan, 
positional differences that are multiples of the grid period 
in the direction of scanning will not produce measurable 
differences in the detector signal. This 'grid-step ambigu- 
ity' is normally resolved by the combination of scans from 
a variety of directions in the course of the mission. 

1 Following the convention of the Hipparcos Catalogue, an 
asterisk (*) is used to denote a true arc, as in p a * = p a cos 8. 
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2.6. Changing the reference point 

The choice of reference point for the Transit Data is ar- 
bitrary, as long as it is sufficiently near the target for 
linearization errors to be negligible f60 arcsec; see 
Sect. 2.5). As mentioned previously, the adopted refer- 
ence point usually corresponds to the H1C values. The 
user may ask why the final astrometric parameters in the 
Hipparcos and Tycho Catalogues were not used instead. 
There were several reasons for this. First of all, the TD 
were derived directly from intermediate files of the NDAC 



reduction (Sect. 2.1), which were compiled in their final 
version almost a year before the astrometric parameters 
of the Hipparcos Catalogue became available. Secondly, 
most of the objects are double or multiple, and it is not 
evident which point to use in such cases. Finally, TD are 
available also for cases where no valid astrometric solution 
is provided in the Hipparcos Catalogue. In those cases it 
would have been necessary to use something like the HIC 
data anyway. Although it would have been possible, using 
the formulae below, to transform the TD to some other, 
possibly more natural or desirable reference point, such a 
process would in the end have been largely arbitrary. In 
addition, since any manipulation of the data always entails 
some risk of error, however small, it was felt better to leave 
the reference points as defined in the NDAC intermediate 
data files. The only modification applied was the transfor- 
mation from the intermediate NDAC reference frame into 
the ICRS system, which was deemed essential; this was 
achieved by rotating the positions and proper motions of 
the reference points, rather than modifying the TD coeffi- 
cients. (This explains why the position and proper motion 
of the reference point never coincide exactly with the HIC 
values.) 

In some circumstances it may be desirable to change 
the reference point for a set of TD. This is particularly 
relevant for the construction of aperture synthesis images 
(Sect. where the reference point defines the centre 
of the image. An object with a poorly determined pre- 
Hipparcos position may be severely offset from the centre 
of the reconstructed map, possibly falling outside the map 
altogether or creating a false image due to aliasing from a 
position outside the map. Not only the positional offset, 
but also errors in the proper motion and parallax of the 
reference point may create problems for the image recon- 
struction. The effect of such errors will be a blurring of 
the image due to the relative motion between the refer- 
ence point and the true object. The proper motion error 
produces a blurring along a straight line in the map while 
the parallax error causes an additional elliptical blurring. 
Most objects with large proper motions (> 100 mas yr" 1 ) 
or parallaxes (> 100 mas) have reasonable, non-zero esti- 
mates of these quantities in the Input Catalogue, so the 
blurring effect is usually not very serious for qualitative 
evaluation of the images. However, for a quantitative anal- 
ysis the effect needs to be considered. 



Changing the astrometric parameters of the reference 
point requires that the phase of each transit is adjusted 
to take into account the apparent positional offset of the 
new reference point from the old one at the time of the 
transit. Fortunately, this is easily done by means of the 
spatial frequency components f x , f y , f p provided for each 
transit. Let Aa* = (a' — ao)cosao, A<5 = 5' — <5o , etc., 
be the (small) differences between the new (') and the old 
reference point in terms of the five astrometric parameters. 
The required phase shift for a particular transit is then 
given by Eq. (||) . The corresponding transformation of the 
Fourier coefficients in Eq. (j|) is 

b\ = h , 

b' 2 = 62 cos 4> — 63 sin <j) , 

b' 3 — 62 sin 4> + &3 cos <p , 

b' 4 = 64 cos 2<f> — 65 sin 20 , 

b 5 = 64 sin 2<fi + 65 cos 20 . 



(7) 



An example of the change of reference point is shown in 
Figs. I and @. 

2. 7. Multiple pointings and target positions 

One complication in the Hipparcos satellite operation and 
data reductions was caused by double and multiple stars 
having separations roughly in the range 10 to 30 arcsec. 
As already mentioned, the main Hipparcos detector had 
a sensitive area (IFOV) of about 30 arcsec diameter. This 
area could be directed towards any pre-defined point on 
the sky currently within the telescope field of view. The 
IFOV pointing, calculated from the real-time knowledge of 
the satellite attitude and the celestial position of the tar- 
get, typically had errors of 1 to 2 arcsec rms. Ideally, no 
star should be observed while its image was close to the 
edge of the IFOV, where the guiding errors might pro- 
duce a distorted signal. For double and multiple systems 
with separations less than about 10 arcsec this could be 
achieved by centering the IFOV somewhere in the mid- 
dle of the system, so that all components remained within 
the flat-topped central part of the IFOV sensitivity profile. 
For systems with separations greater than about 30 arc- 
sec the individual components (or subsystems with sepa- 
rations below 10 arcsec) could be observed as single stars, 
again avoiding signal distortion from the IFOV edges. 

However, systems with intermediate separations (~ 10 
to 30 arcsec) could not be observed without some adverse 
effects of the IFOV edges. In order to allow at least some 
useful astrometric information to be extracted for such 
systems, each component (or subsystem) received a sep- 
arate pointing. For example, HIP 70 and HIP 71 formed 
such a two-pointing system with a separation of about 
15 arcsec. When pointing at the brighter star (HIP 71, 
Hp ~ 8.4), the other component (HIP 70, Hp ~ 10.6) 
would be just at the edge of the IFOV and a (variable) 
fraction of its signal was added to that of the brighter 
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star. Conversely, when pointing at the fainter component, 
some fraction of the brighter star's signal would be added. 
Proper reduction of such systems must consider the mu- 
tual (and possibly distorted) influence of each component 
upon the other, as was indeed done in the Hipparcos data 
reductions. In some multiple systems three different point- 
ings were needed. 

The situation is further complicated by the fact that 
the targeted position of the IFOV was sometimes updated 
in the course of the mission, usually because the original 
(ground-based) position was found to be wrong by several 
arcsec. Knowledge of both the original and the updated 
position, and the time of updating, may then be necessary 
for proper interpretation of the observations. 

In order to cope with two- and three-pointing systems 
as well as updated positions, the concept of 'target posi- 
tions' was introduced in the TD. For a set of TD referring 
to a particular double or multiple system, target positions 
are defined by their offset coordinates (Ace*, AS) from the 
adopted reference point, rounded to the nearest arcsec. 
In most cases there is just one target position, coincid- 
ing with the reference point [offset coordinates = (0,0)]. 
For multiple-pointing systems there is at least one target 
position for each pointing, with a different HIP number 
attached to each pointing. For objects whose coordinates 
relative to the reference point changed in the course of the 
mission, a new target position was introduced whenever 
the offset coordinates changed by more than 1 arcsec. To 
within the errors of the real-time attitude determination 
(normally 1-2 arcsec rms) it can therefore be assumed that 
the IFOV was pointed to the specified target position. 

In the TD the results from different pointings per- 
taining to the same system have been collected together 
and expressed relative to a common reference point. This 
has been done for all systems deemed to be 'difficult' in 
the sense explained above, but not for very wide systems 
where the mutual influence of the component signals was 
negligible. The TD moreover contains the bookkeeping 
data necessary to calculate the actual target positions in 
each transit. 



3. The Transit Data files 

This section gives a rather detailed description of the con- 
tents of the TD files, complementing the formal descrip- 
tion in Vol. 1 of ESA ( |1997| ) and providing additional 
explanation of the data items. All relevant data are con- 
tained in two ASCII files, both located on Disk 6 in Vol. 17 
of ESA ( p97p : the TD index file (hip.j.idx, ~ 1 Mb) and 
the TD file (hip_j.dat, ~ 553 Mb). A summary of the 



Table 1. Summary of the contents of the TD index and 
data files on Disk 6 of the Hipparcos Catalogue. Data fields 
in the Transit Data File are separated by the character 
'I'. In the header record there are 13 fields designated JH1 
through JH13, etc. 

Transit Data Index File (hipj.idx) 
120 416 records, each of 7 bytes+CR+LF = 9 bytes 
The nth record in this file points to the header record 
for HIP n in hip_j.dat (Fig. 2) 

Transit Data File (hip_j.dat) 
4 351 156 records, each of 125 bytes+CR+LF = 127 bytes 
The data for a given entry consist of one header record, 
one pointing record, and Nt transit records. 

Header Record: 
JH1-3 HIP identifiers (JH2-3 = if not used) 
JH4 Np = number of target positions 

JH5 Nt = number of transit records 

JH6-10 ao, So, 7ro, p a *o, fJ-so = reference point 
JH11-13 assumed colour indices (V—I) for JH1-3 

Pointing Record: 
JP1 index (1 to 3) to JH1-3 for target position 1 

JP2 offset in a for target position 1 (arcsec) 

JP3 offset in 8 for target position 1 (arcsec) 

JP4-6 same as JP1-3 but for target position 2 

JP25-27 same as JP1-3 but for target position 9 

Transit Record: 

JT1 Ip = target position (1 to Np) for the transit 

JT2 t — epoch of the transit (years from J1991.25) 

JT3-5 fx, f y , f P = spatial frequencies 

JT6 ln&i 

JT7-10 bi/bi, i = 2. . .5 

JT11-15 Incrj, i = 1 . . . 5 (standard errors of bi) 

JT16-17 si, S2 = colour correction factors 

JT18 (7 at t = additional attitude noise (mas) 

JT19 computed (0) or assumed (1) standard errors 



3.1. The index file: hip_j.idx 

The TD index file is included to facilitate accessing the 
TD file. It contains a pointer from the HIP number to the 
corresponding record in the TD file where data on that 
object can be found. Since HIP numbers range from 1 
to 120416, there are exactly 120416 records in the index 
file. Each record consists of a 7-character integer, which is 
the record number in the TD file for the relevant header 



contents of the two files is in Table [l], while important 
relations among the data are illustrated in Fig. [| 



record (Sect. 3.2.1). For example, TD for the double star 
HIP 7 can be accessed by first reading record number 7 in 
the index file (hip_j.idx). The content of that record is the 
integer 129. Record number 129 in the TD file (hip_j.dat) 
is thus the header record for HIP 7, and information on 
that object is contained in that and subsequent records 
of the TD file. If no TD are available for a given HIP 
number, then the corresponding index file record contains 
the number — 1. 
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hipj.idx 

Rec # Index 
(= HIP #) data 



hipj.dat 



Transit Data 
3 1 126 .00497547... 

^ioo o o o o o o... 

1 -1.3271349 696039 -817248... 

1 -1.3270956 694424 -818615... 

1 -1.3268913 693072 -819757... 
ETC. 



420 
421 
422 
423 
424 
425 



129 
248 
-1 



17238 
-1 
17057 
17238 
17581 




1 1.8786947 -376914-1005176.. 



7 1 117 .02236934... 
■ ]! 0... 
- 1 -1.2848122 -376914 -10051.. 
ETC. 



17238 
17239 
17240 
17241 
17242 
17243 
17244 







424 421 


2 341 1.29074007.. 



■\2 -18 -12 1 0... 



1 -1.3237247 637663 -863563... 

2 -1.3237246 637667 -863559... 
2 -1.3236856 636772 -864213... 
1 -1.3236855 636772 -864216... 
1 -1.3234815 634086 -866191... 

ETC. 



HIP 3 header 
Pointing record 
Transit record 1 
Transit record 2 
Transit record 3 



Transit record 126 
HIP 7 header 
Pointing record 
Transit record 1 



HIP 421/424 header 
Pointing record 
Transit record 1 
Transit record 2 
Transit record 3 
Transit record 4 
Transit record 5 



Fig. 2. Relations among the contents of the index and 
data files. 



For two- and three-pointing systems (Sect. 2.7) there 
are two or three different index file entries pointing to the 
same TD record. An example is shown in Fig. |^, where 
the index file entries for HIP 421 and HIP 424 both point 
to the 17238th record in the TD file. 

For efficient accessing of the TD it is recommended 
that the index file is read into computer prime memory as 
a one-dimensional integer array of length 120416. This ar- 
ray is then used as a look-up table for finding data in 
the TD file. The TD file has a constant record length 
of 127 bytes (including the two end-of-record characters 
\r\n = CR+LF), which allows simple direct access to any 
given record number. 

3.2. The TD file: hip_j.dat 
3.2.1. The header record 

The header record is the first record in the TD file contain- 
ing data on a specific system. The index file always points 
to this header record for a given HIP number. The header 
record contains general information about the system in 
13 data fields. In holding to the Hipparcos Catalogue con- 
ventions, the 13 data fields in the header record will be 
called JH1, JH2, . . . , JH13. 



The first three data fields (JH1-3) contain the HIP 
numbers relevant to the subsequent transit records. JH1 is 
always defined; JH2 and JH3 are only defined for two- and 



three-pointing systems (Sect. 2.7). The number of point- 
ings, or more accurately the number of different target po- 
sitions (Np), is specified in field JH4. Details on the rela- 
tive pointings are given in the pointing record (Sect. |3T2~2 ). 

JH5 gives the number of transits (Nt), or observations, 
for a given system. Since there is one record for each tran- 
sit, this is also the number of transit records. The total 
number of records for an object is then Nt + 2, where 
the header and pointing records make up the extra two 
records. 

JH6-JH10 give, in order, the right ascension (deg), 
declination (deg), parallax (mas), proper motion in right 
ascension (mas yr" 1 ) and proper motion in declination 
(mas yr _1 ) of the reference point. The Fourier coefficients 
&1-&5 of the subsequent transit records are expressed rela- 
tive to this reference point. It is important to realize that 
these values are usually derived from the Input Catalogue 
and therefore do not agree with the astrometric parame- 
ters given in the Hipparcos Catalogue. For instance, the 
parallax of the reference point is often zero. The posi- 
tion refers to the epoch J1991.25, which was close to mid 
mission. The reference system is the same as that of the 
Hipparcos Catalogue, viz. the Internationa l Cele stial Ref- 
erence System (ICRS) (Feissel & Mignard |l99^) . 

JH11-JH13 hold the assumed colour index V— I (mag) 
for each of the three HIP numbers in JH1 through JH3. 
This value V — I was used in rectifying the signals from 
each transit. 

3.2.2. The pointing record 

The pointing record defines the offset of the target position 
of each subsequent transit with respect to the reference 
point. The pointing record may define several different 
target positions as required for multiple-pointing objects. 
The actual number of different target positions (Np) is 
given in JH4 of the header record. Each target position 
is specified by three numbers: the first takes the value 1, 
2 or 3, depending on which of the HIP numbers in JH1, 
JH2 or JH3 that the target position refers to. The second 
and third numbers give the offset in a and S, respectively, 
of the target position from the reference point. These last 
two values are rounded to the nearest arcsecond. 

Because every record of the TD file has 125 bytes, this 
allowed up to nine different target positions to be defined 
in the pointing record. In reality the maximum Np was 7. 

Most objects used a single pointing, and the astromet- 
ric values in the Input Catalogue were sufficiently accurate 
that no updating was required. Moreover, the Input Cat- 
alogue values were usually taken as the reference point for 
the transit data. In this case Np = 1 and the first and only 
entry in the pointing record contains the three values: '1 
0' (cf. the excerpts for HIP 3 and 7 in Fig. |). The '1' 
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refers to the HIP number in JH1. The '0 0' means that 
the target position coincided (to the nearest arcsec) with 
the reference point. Subsequent transit records have '1' in 
the first field (JT1) indicating that the target position for 
each transit was as given in the first entry of the pointing 
record. 

In the trivial case of a single, well-centered pointing, 
as illustrated by HIP 3 and 7 above, the pointing record 
is really not needed. However, for those objects with mul- 
tiple HIP entries and pointings, the pointing data provide 
important information. This is exemplified by the excerpt 
for the two-pointing system HIP 421+424 shown in the 
lower part of Fig. |2[ As indicated in the header record, 
there are 341 transits for this system, made using two dif- 
ferent target positions. The first transit for this system 
(in record 17240) was made with the IFOV pointing at 
the first target position ('2 -18 -12'), i.e. referring to 
the second HIP entry (421) and pointing 18 arcsec to the 
west and 12 arcsec to the south of the reference point. The 
next transit (record 17241) was made in the second point- 
ing ('1 0'), i.e. referring to the first HIP entry (424) and 
centred on the reference point. The third transit (record 
17242) was also made in the second pointing, while the 
fourth transit was made in the first pointing, and so on. 

The fact that Hipparcos had two fields of view, sepa- 
rated by the 'basic angle' of 58°, is largely irrelevant for 
the use of the TD. There is in fact no flag in the TD telling 
in which field of view a particular transit occurred. 

3.2.3. The transit data record 

The actual scan information is contained in the transit 
data record. This section will cover the entries contained 
in these records, as well as some physical interpretations. 

The first field (JT1) for each record tells which of the 
Np (= JH4) target positions was used to observe the tran- 
sit. The second field (JT2) gives the time of the observa- 
tion, in years, from the epoch J1991.25. The time is mea- 
sured in Julian years of exactly 365.25 days. The epoch 
J1991.25 is equivalent to the Julian date JD 2 448 349.0625 
on the Terrestrial Time (TT) scale. 

The next three entries (JT3 to JT5) contain the spa- 
tial frequencies f x , f y and f p defined as in Sects 2.4 and 
2.5. They are expressed in units of rad rad" 1 (radian of 



modulation phase on the grid per radian on the sky). 

The entries JT6 through JT10 contain, in coded form, 
the five Fourier coefficients 61-65 defined by Eq. (1). JT6 
contains the natural logarithm of b\ (In 61) and JT7 to 
JT10 contain the normalized coefficients (62/61, 63/61, 
64/61 and 65/61). This coding was adopted because 61, 
being the mean intensity of the signal, is always positive 
and spanning a wide range of values according to the mag- 
nitude of the object; 62 through 65, on the other hand, may 
have either sign but are always numerically smaller than 
61. 



The next five entries (JT11-15) contain the natural 
logarithms of the standard errors (o\-<j§) for each of the 
Fourier coefficients. JT19, the last entry in this record, 
is a flag for the standard errors. JT19 is normally set to 
zero. It is set to '1' if the standard errors could not be 
computed in the normal way, but were estimated roughly 
from the standard errors of the other transits. A non-zero 
value in JT19 is thus a warning that the standard errors 
in JT11-15 should be treated with caution. 

JT16-17 contain two colour correction factors, s\ and 
S2- As previously mentioned, the TD are rectified signals, 
i.e. corrected for calibrated variations in the geometric and 
photometric responses. Among the many factors that have 
been taken into account in this process are the colour ef- 
fects with respect to position in the field of view and time 
variations of the photometric calibrations. These correc- 
tions depend on an assumed colour for the object, as given 
in JH11-13. If it should turn out that the assumed colour 
for a particular transit was considerably in error, si and 
S2 provide a possibility to correct (approximately) the 
Fourier coefficients accordingly. See Sect. 2.9.3 in Vol. 1 
of the Hipparcos and Tycho Catalogues for details of the 
correction procedure. 

The second to the last entry (JT18) deals with the 
errors in the attitude determination. The phase determi- 
nation of each scan was basically limited by photon noise, 
but the standard errors in JT11-15 also contain a contri- 
bution from the along-scan attitude uncertainty. It turned 
out to be very difficult to treat the attitude errors rigor- 
ously in the TD, and as a result these errors were gen- 
erally underestimated. JT18 contains an estimate of the 
additional attitude noise (cr a tt, in mas) required in order 
to derive, from the TD, astrometric standard errors that 
are consistent with the main processing of Hipparcos data. 

3.2.4. Known errors in the Transit Data 

After publication of the Transit Data CD-ROM, six cases 
of Fortran format overflow errors have been discovered in 
the file hip_j.dat. No useful data are lost because of the 
errors, but special measures may be needed to read the 
corresponding records. The interface programs described 
in Sect. [| automatically correct these errors. Further de- 
tails can be found at the Internet address given in Sect. § 

4. Aperture synthesis imaging 

The Hipparcos satellite was not designed for imaging and 
did not contain any imaging device such as a CCD camera. 
The combination of a modulating grid and the IDT, while 
well adapted to the observation of isolated point sources, 
was far from ideal for the observation of more complex 
resolved objects. As explained in Sect. [2^ the grid es- 
sentially extracted two spatial frequencies (with periods 
1.2074 arcsec and 0.6037 arcsec in the direction of each 
scan) of whatever intensity distribution was within the 
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30 arcsec sensitive spot. This is much more reminiscent of 
spatial intcrfcromctry than of normal optical imaging. In 
fact, the two harmonics of the detector signal correspond 
to the fringes produced by an object in two interferom- 
eters with baselines of 94 mm and 188 mm, respectively 
(assuming an effective wavelength of 550 nm). Image re- 
construction techniques using intcrferometric observations 
have for a long time been standard in the radio astronom- 
ical community. It was therefore quite natural to apply 
these techni ques to the Hipparcos data (Lindegren 1982 
Quist et al 



1997|) 



4-1- Mathematical foundation 

Equation ([!]) gives the signal for a point source located at 
the position r = (x,y) relative the reference point. Now 
consider an extended object with the general brightness 
distribution B{r). Summing up the contributions to the 
detector signal from each element of the sky we find (apart 
from a constant scaling factor) 



J(p) = J J S(r)B(r)[l + Mi cos(p + f ■ r) 
+ M 2 cos(2p + 2/ • r)] d 2 r 



(8) 



where S(r) is the sensitivity profile of the instantaneous 
field of view. Introducing the Fourier transform 



V(f)= // S(r)B(r)exp(i/T)d 2 r 



(9) 



we find that Eq. (j|) can be written 

I(p) = 7(0) + Re[V(f)]Mi cosp - lm[V{f)]Mi sinp 

+ Re[U(2/)]Af 2 cos2p- im[F(2/)]M 2 sin2p . (10) 

Comparison with Eq. (Q) shows that 

V(0) = h , 

V(f) = (b 2 - i&aVMi , V(-f) = V*(f) , 
V(2f) = (b 4 -ib 5 )/M 2 , V(-2f) = V*(2f), (11) 

where the asterisk denotes the complex conjug ate. Since 
M\ and M 2 are conventional constants (Sect. |2.3| ) it is 
seen that a single transit defines the complex function V 
in the five points 0, ±/, ±2/, of the spatial frequency 
plane. However, the conjugate symmetry of V means that 
there are only three independent complex visibilities per 
transit. Successive transits of the same object are made 
at different spatial frequencies / = {f x ,fy) and in the 
course of the mission knowledge of the function V is built 
up in a number of different points. To the extent that 
S(r)B(r) remains constant over the mission, it may then 
be recoverable from V using standard image reconstruc- 
tion techniques. 

In the context of radio interferometry and aperture 
synthesis, we may identify V(f) with the complex visi- 
bility function associated with the source brightness dis- 
tribution B{r) and single-antenna reception pattern Sir) 



(Thompson et al. 1994). The visibility function is usually 
expressed in terms of coordinates (it, v) which give the 
projection of the interferometer baseline on the sky plane 
and are expressed in wavelengths. The relation to the TD 
spatial frequency components is simply 



/»/27T, 



f y /2w. 



(12) 



Our reference point is equivalent to the phase reference 
position used in connected-element radio interferometry 
(Thompson et al. 1994), or to the strong reference point 
source (typically a quasar) used in phase-referenced VLBI 
observations (Lestrade et al. 1990). 

The distribution of the observations in the uv plane 
is all-important for the possibility to reconstruct compli- 
cated images from the measured visibilities V(u, v). Un- 
fortunately the Hipparcos scanning law and the use of a 
modulating grid with just a single period seriously limit 
the uv coverage of the TD. According to Eqs. (||) and fll2| ) 
the coverage is limited to the central point (u, v) = (0, 0) 
and two concentric rings with radii ~ 170 830 and 341 660 
wavelengths. Moreover, for objects in the ecliptic region of 
the sky (ecliptic latitude \/3\ Si 45°) the scanning law con- 
strains the scan angle 9 such as to produce a gap of 'miss- 
ing' scans roughly in the east-west direction. At \(3\ ~ 47° 
there is instead a surplus of scans in the east-west di- 
rection (cf. Fig. ||). For high-latitude objects, finally, the 
coverage is usually more uniform in 9. 

In continuing the analogy with radio interferometry, 
we will discuss in the next section how images can be 
produced using the Transit Data. Utilizing the experi- 
ence developed for aperture synthesis imaging, we use only 
publicly available software for producing and deconvolving 
these images. 

4-2. Conversion of TD to UV- FITS format 

The images presented here were produced using the Cal- 
tech Difmap software package (Shepherd 1997 ). Another 
software package that could be used for the analysis is 
AIPS (Astronomical Image Processing System) developed 
at the National Radio Astronomy Observatories (NRAO). 
Both Difmap and AIPS take input data in the form of 
FITS files, using the 'Random Groups' format (NOST 



1994). This amendment to the basic FITS format is more 



commonly called the UV-FITS format, since it is used al- 
most solely for interferometry data. 

We have written and made publicly available (Sect. |6|) 
a concise Fortran program which reads TD from the CD- 
ROM format (hip_j.idx and hip_j.dat) for any given object 
and produces an output file in the UV-FITS format. In the 
following we describe some of the features of the UV-FITS 
format and how the TD were adapted to it. 

The basic FITS (Flexible Image Transport System; 
Wells et al. [1981[ NOST |1993j ), well known in optical as- 
tronomy, was designed to transport digital data in the 
form of rt-dimensional regular arrays, such as CCD images, 
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with associated information on coordinates, dates, scales, 
units, etc. given in an ASCII header. However, aperture 
synthesis visibility data do not come in regular arrays, at 
least not in all axes, and thus an amendment to the ba- 
sic FITS was required allowing the definition of ordered 



sets of small arrays (Greisen & Harten 1981). In the fol- 
lowing we assume that the reader has some rudimentary 
knowledge about the basic FITS format. 

In a UV-FITS file each random group contains the 
visibility data associated with a particular point in uv 
space and time. The group consists of a set of parame- 
ters followed by a regular array of measurements. The pa- 
rameters are, for instance, the uv coordinates and date of 
the measurements. The measurement array may be multi- 
dimensional with, for instance, the different frequencies 
and Stokes (polarization) components marked along two 
of the axes. In the UV-FITS header the use of random 
groups is signified by having a first axis of length zero 
(NAXIS1 = 0) and by setting the keyword GROUPS = T 
(true). Visibility data are stored as three values in the 
measurement array, namely the real part of the visibility, 
the imaginary part, and an associated weight. In the array, 
this corresponds to an axis of length 3 and type COMPLEX. 

For the Hipparcos TD there are three visibility mea- 
surements per transit, corresponding to the spatial fre- 
quencies 0, / and 2f. The total number of random groups 
is therefore GCOUNT = 3V T . PCOUNT = 6 parameters spec- 
ify each group, namely the Fourier coordinates (u,v,w), 
the baseline (defining which pair of antennae that formed 
the interferometer), and the date of the observation (split 
in two numbers containing the integer and fractional parts 
of the Julian date). For the TD the w coordinate is always 
zero. 

The UV-FITS format requires that the (u, v) coordi- 
nates are expressed in seconds, while in Eq. ( |l2[ ) they are 
dimensionless. The conversion factor requires the specifi- 
cation of a reference wavelength, for which we arbitrar- 
ily adopted Ao = 550 nm. The corresponding scale factor 
for the (u, v) coordinates in Eq. (jlj) is then PSCAL1 = 
PSCAL2 = Ao/c, where c is the speed of light. For consis- 
tency, the frequency associated with each TD observation 
must then be given as c/Ao. 

The measurement array in each group is 5-dimensional 
(NAXIS = 6, since the first axis has zero length for group 
data). Its size is NAXIS2 x NAXIS3 x NAXIS4 x NAXIS5 x 
NAXIS6 = 3x1x1x1x1 = 3, where the axes are of 
type COMPLEX, STOKES, FREQ, RA and DEC, respectively. The 
values on each axis are specified in the header by CRVALn 
for n — 2 ... 6; in particular the frequency is given by 
CRVAL4 = c/Ao and the reference position by CRVAL5 = ao 
and CRVAL6 = 5q. The three values in the measurement 
array are, as mentioned before, the real and imaginary 
parts of V(u, v) and an associated weight. For the TD the 
weight is always equal to 1. 

The aperture synthesis programs also require the 
names and geocentric positions of the antenna stations 



/' 


1 1 


' 1 1 


' / 

I 


/' 

/ 

i 








/ 1 

1 1 

! 

t 









0.2 





U (10 6 \) 



-0.2 



Fig. 3. UV coverage for HIP 97237. For every point (u, v) 
marked in this diagram there is a measurement of the com- 
plex visibility V representing the amplitude and phase of 
the detector signal during a certain transit. Each transit 
provides five equidistant visibilities (including the origin) 
oriented along a straight line through the origin. The ce- 
lestial orientation of the line is normal to the slits of the 
Hipparcos modulation grid at the transit, u is the spatial 
frequency in the east-west direction, expressed in modu- 
lation periods ('wavelengths') per radian on the sky; v is 
the north-south component of the spatial frequency in the 
same unit. 



to be specified, although this is rather pointless in our 
case. We formally specify six stations and identify a dif- 
ferent pair with each spatial frequency. The stations are 
arbitrarily named to form the acronym 'HIPUVF', which 
will appear on some plots. 

An apparent limitation of the FITS format is the lack 
of keywords for proper motion and parallax. Until such 
keywords become standard, the proper motions and par- 
allax values of the reference point are given in the ASCII 
header as comments. 



4-. 3. Example images 

Once the data are converted into a suitable form, im- 
ages are quickly produced using standard aperture syn- 
thesis programs. Using the Difmap package from Caltech 
(Sect. ^), we give here as an example a reduction of the TD 
for HIP 97237. In the Hipparcos Catalogue, no astrometric 
solution is given for this relatively faint (V ~ 12.4) object. 
In the Hipparcos Input Catalogue it is noted as a double 
star (CCDM 19458+2707) with separation 0.9 arcsec and 
component magnitudes 12.7 and 13.6. The ecliptic latitude 
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Fig. 4. The dirty beam (point spread function) associated 
with the UV coverage of HIP 97237 shown in Fig. |. 

Fig. 5. The dirty map of HIP 97237, obtained by an in- 
verse Fourier transform of all the measured complex visi- 
bilities. 

of the object is (3 ~ +47°. The reference point for this ob- 
ject is a = 296.43849356 deg, S = +27.12735771 deg, 
7To = 97.00 mas, /i Q *o = —25.43 mas yr _1 , figo = 
-1228.86 mas yr _1 . 

Figure | shows the UV-coverage of HIP 97237. It is 
seen that the object received rather many scans, although 
predominantly in the (ecliptic) east-west direction. The 
'dirty beam' (Fig. ||) reflects this anisotropy as a charac- 
teristic pattern along the east-west section, reminiscent of 
the basic light modulation curve in Fig. |[ 

The 'dirty map', basically obtained as the inverse 
Fourier transform of the complex visibilities, is shown in 
Fig. |^. Already from this image it is obvious that the 
object was offset from the expected (HIC) position by 
about 5 arcsec. This is probably the reason why an ac- 
ceptable solution was not found in the Hipparcos astro- 
metric reductions. Deconvolution of the dirty image, us- 
ing the dirty beam as kernel (point spread function), can 
be achieved by means of the CLEAN algorithm (Hogbom 



1974) implemented in Difmap. One result of this pro- 



cess (which depends on several parameters selectable in 
Difmap) is shown in Fig. Ifj. Both components of the dou- 
ble star are now clearly seen. The offset of the primary 
component from the reference point, estimated from the 
cleaned image, is (Aa*,A5) ~ (+3.6,-1-3.6) arcsec. The 
position of the secondary component relative to the pri- 
mary is approximately 0.9 arcsec towards position angle 
330° . The power in the primary peak of the cleaned map 
is 0.0455 units, where 6200 units corresponds to magni- 
tude Hp — (Sect. (2~3| ); the magnitude of the primary 
can thus be estimated at Hp ~ 12.8. 

Figure [?] illustrates the cha nge of reference point for 
the TD described in Sect. |2.6| . The astrometric param- 
eters of the primary component in HIP 97237, relative 
to the reference point, were estimated by means of the 
model fitting procedure described in Sect. |[ The ap- 
proximate results were (cf. Fig. ||) Act* = +3603 mas, 
AS = +3640 mas, An = — 16 mas, A/i Q * = —50 mas yr _1 , 
A us = — 25 mas yr _1 . Applying the corresponding phase 
shifts to the TD, according to Eqs. (0) and (||), effectively 
changes the reference point to coincide with the primary 
component. Performing the image synthesis on the mod- 
ified TD gives the cleaned image in Fig. [j]. As expected, 
the primary now appears at the centre of the map. The 
shift in parallax and proper motion of the reference point 
improves the relative phasing of the superposed scans, re- 
sulting in a slightly increased peak power (from 0.0455 to 
0.0457 units). 



Fig. 6. The cleaned map of HIP 97237, obtained by de- 
convolution of Fig. H with the beam in Fig. [|. The cleaned 
map clearly reveals the two components of the double star. 

Fig. 7. The cleaned map of HIP 97237 resulting from mod- 
ified TD, in which the astrometric parameters of the ref- 
erence point were shifted to coincide with the astrometric 
parameters of the primary component, determined using 
the model fitting procedure described in Sect. ^. 

5. Model fitting 

The aperture synthesis imaging attempts to reconstruct 
the brightness distribution on the sky, in principle with- 
out making any a priori assumption about the object. This 
is excellent for exploring cases where the nature of the ob- 
ject is uncertain, e.g. concerning the number of resolved 
components in a multiple star or their approximate posi- 
tions. The method is less useful for accurate quantitative 
evaluation, in particular because the offsets in parallax 
and proper motion merely produce a blurring of the im- 
age. Since the objects of interest here consist of a small 
number of point sources, direct modeling of the TD in 
terms of simple superposed signal components is usually 
possible. Such model fitting provides the most direct and 
accurate estimates of specific object parameters such as 
the trigonometric parallax or orbital elements. In this sec- 
tion we outline the fitting procedure and give an example 
of its practical realization by means of a publicly available 
computer program. 

5.1. General method 

Let us assume that the object consists of n point sources 
with intensities Aj and positions Xj , yj relative the refer- 
ence point (j = 1 . . . n) . In general Aj , Xj , yj vary with 
time, and so may be different for the different transits 
of the same object. Given a specific model of the object 
we express Aj, Xj, yj as functions of time t and a set of 
model parameters a. For instance, in the case of a non- 
variable orbital binary, a would consist of 15 parameters, 
viz. the five astrometric parameters of the mass centre, 
the magnitude of each component, the mass ratio, and 
seven elements for the relative orbit. Generally speaking, 
the object model is thus completely specified by n and 
the functions Aj(t, a), Xj(t, a), yj(t, a) for j = 1 . . . n. In 
the equations below we suppress, for brevity, the explicit 
dependence on t and a. 

For a given transit the expected signal is modeled as 
the sum of the signals from the individual components, 
using Eqs. (|) and (§). Thus, 



Ml CO&(pk + fxXj + fyVj) 

M 2 cos 2(p k + f x Xj + f y yj)} . 



(13) 
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Expanding the trigonometric functions and equating the 
terms with those in Eq. (Q) yields 

h = J2 A i> 

j 

b 2 = Mi ^2 A.j cos(f x Xj + fyyj) , 

3 

6 3 = -M 1 ^2A d wa.{f x x j +f v y j ), (14) 

3 

j cos2(f x Xj + fyUj) , 

j 

b 5 = ~M 2 ^ A 3 Sin 2 ifxXj + f y yj) ■ 



Recall that Aj, Xj, yj depend on the model parameters a. 
The general procedure is then to adjust a in such a way 
that, for the whole set of transits, the calculated signal 
parameters 61-65 from Eq. ( |l4| ) agree, as well as possi- 
ble, with the observed values. The adjustment may use 
the weighted least-squares method, using the standard er- 
rors of the observed signal parameters to set the weights; 
but other (and more robust) metrics can also be used. In 
general the problem can be formulated as a constrained 
minimization problem in the multi-dimensional model pa- 
rameter space. 

The trigonometric functions in Eq. ( |l4| ) mean that the 
signal parameters bi depend in a highly non-linear man- 
ner on the model parameters which affect Xj and yj. For 
instance, in terms of a displacement of one of the point 
sources, the effect on 64 and 65 is approximately linear 
only for displacements less than about 1/2/ ~ 0.1 arcsec, 
corresponding to 1 rad change in the modulation phase. In 
the aperture synthesis imaging this non-linearity is mani- 
fest in the complex structure of the 'dirty beam' (Fig. [|) at 
all spatial scales larger than about 0.1 arcsec. Additional 
non-linearities in the complete object model may result 
from the geometrical description of the source positions, 
e.g. in terms of orbital elements. 

The non-linearity of the object model has two impor- 
tant consequences for the model fitting. Firstly, it is usu- 
ally necessary to use a non-linear, iterative adjustment al- 
gorithm, such as the Levenberg-Marquardt method (Press 
et al. 1992 ). Secondly, a good initial guess of the model pa- 
rameters is usually required. In particular the parameters 
directly affecting the positions of the point sources need to 
be specified to within (what corresponds to) a few tenths 
of an arcsec. Without a good initial guess, the adjustment 
algorithm is likely to converge on some local minimum, 
typically resulting in positional errors of (approximately) 
an integer number of grid periods. The correct solution, 
corresponding to the global minimum, may in principle al- 
ways be found through sufficiently extensive searching of 
the parameter space. Alternatively, sufficiently good ini- 
tial guesses of the point source positions can often be ob- 
tained from the aperture synthesis imaging. 



Various least-squares model fitting procedures were 
used for the reduction of double and multiple stars dur- 
ing the construction of the Hipparcos Catalogue (see 
Mignard et al. 1992 and references therein). The double- 
star processing of the NDAC data reduction consortium 
(Soderhjelm et al. 1992| ) essentially used the technique 
outlined above, taking the so-called Case History Files (a 
precursor to the TD) as input. 

Perhaps the greatest potential of the TD lies in the 
possibility to combine the Hipparcos data with indepen- 
dent observations from other instruments and epochs. For 
instance, full determination of a binary orbit generally 
requires data covering at least a whole period. Ground- 
based speckle observations can sometimes provide this, 
constraining the geometry of the relative orbit much bet- 
ter than the Hipparcos data alone, and in turn leading 
to a better-determined space parallax. In some favourable 
cases the location of the mass centre in the relative orbit 
(and hence the mass ratio) can be determined (Soderhjelm 
et al. |1997| ; Soderhjelm |1999| ). 

One complication of the Hipparcos double star process- 
ing has been the wide variety of applicable object models, 
and the consequent need to experiment and interact with 
the solutions. This process may be much facilitated by 
using general and flexible software for the model fitting, 
rather than highly specialized routines. An example of this 
is given below. 

5.2. Model fitting using GaussFit 



GaussFit (Jefferys et al. 1988a, 1988b) is a general pro- 
gram for the solution of least squares and robust estima- 
tion problems, developed as a platform to facilitate astro- 
metric reduction of data from the Hubble Space Telescope. 
It is written in the C programming language and may thus 
be run under a variety of operating systems. In this sec- 
tion we outline the use of GaussFit for model fitting to 
the TD, again using the binary HIP 97237 as illustration. 



GaussFit was used by Soderhjelm ( 1999 ) in a sys- 
tematic re-examination of the solutions for several dou- 
ble and multiple objects, through a combination of TD 
with ground-based observations. Although not illustrated 
in the example below, the introduction of additional data 
(e.g. relative positions from speckle observations) is quite 
straightforward by means of GaussFit. 

To run GaussFit, the user must supply several input 
files. During execution these files are read (and sometimes 
modified) by GaussFit, and additional output files gener- 
ated. For application to the TD model fitting the following 
input files are required. 

— The data file: this contains the observational data, in 
our case the TD. A special program (td2uv.f) is avail- 
able (Sect. ^J) to extract the TD for a given HIP number 
and format them as required by GaussFit. The result- 
ing data file consists of 16 columns and one data line 
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/* mul2_fix.mod - binary model with fixed relative position */ 
observation bl, b2, b3, b4, b5; 

parameter xl, yl, par, xdot, ydot, Al, x2, y2, A2; 
data obs, it, t, fx, fy, fp; 

variable mas = 4 . 8481369E-9, Ml = 0.71, M2 - 0.2485, pi, p2; 

main ( ) 

{ 

while ( import ( ) ) { 

export (bl - (A1+A2)); 

> 

while ( import ( ) ) { 

pi m (fx* (xl+t*xdot) 
p2 m (fx* (x2+t*xdot) 
export (b2 
export (b3 
export (b4 



fy* (yl+t*ydot) 
fy* (y2+t*ydot) 
Ml* (Al*cos (pi) +A2*cos (p2) ) ) ; 
Ml* (-Al*sin(pl)-A2*sin(p2) ) ) ; 
M2* (Al*cos (2*pl) +A2*cos (2*p2) ) ) ; 
export (b5 - M2* (-Al*sin (2*pl) -A2*sin (2*p2) ) ) ; 



fp*par) *mas; 
fp*par) *mas; 



) 



Fig. 8. An example of a double-star model defined in the 
GaussFit programming language. 



per transit. The columns contain a sequential number 
for the transit, the target position index (JT1), the 
time of the transit, the spatial frequencies f x , f y , f p , 
the signal parameters 61-65 and their variances. The 
header of the data file defines the name of the variable 
associated with each column. 

— The model file (cf. Fig. |J): this is a mathematical de- 
scription of the object model written in the GaussFit 
programming language. This language is modeled on 
C, but includes some specific constructs. For instance, 
the declaration of variables distinguishes between 'ob- 
servations' (input data with random errors that need 
to be taken into account in the fitting), 'data' (error- 
free input data), 'parameters' (to be adjusted by the 
program), and ordinary 'variables'. The special func- 
tion import () reads one line of data from the data file. 
The function export (x) sends the equation of condi- 
tion x = to the estimation algorithm, taking into ac- 
count the uncertainties of the observational data that 
went into calculating x. 

— The parameter file: this contains the initial guesses of 
all the model parameters to be estimated. On output it 
contains the estimated parameter values and estimated 
errors. 

— The environment file contains general information 
needed for the model fitting, such as the names of the 
data, parameter and output files; the type of estima- 
tion algorithm to be used (standard least squares or a 
robust method), and stopping rules for the iterations. 

The reader is refer red to the GaussFit User's Manual (Jef- 
ferys et al. 1988b| ) for detailed information. 

Figure || is an example of a GaussFit model file. It 
describes a binary with a fixed positional offset between 
the components (i.e. a long-period binary). The model 
parameters are thus the astrometric parameters of the 
primary relative to the reference point (xl = Aa*i, 
yl = A5i, par = Air, xdot = A/x a *, ydot = A fig), the 



Al = 0.043840289782164 
A2 = 0.0225072066508319 
xl = 3602.736776874136 

xdot = -49.5105454268140330 
yl = 3640.476508541818 

ydot = -24.9326404672069410 
x2 = 4022.0257119596145 
par = -16.2567178488562620 
y2 = 4347.9251385027246 

fit =0.0000034594031099 

End of Iteration 9 . . . 

Sigma Values 



delta 
delta 
delta 
delta 
delta 
delta 
delta 
delta 
delta 

tol = 



0000000022754718 
0000000024041859 
0000111835208815 
0000089562943078 
0000057252486165 
0000142204008415 
0000329749169950 
0000084409108263 
0000428443374577 

0.00001 



Sigma 


Al 


4 


108705745892607e-04 


Sigma 


A2 


4 


11420942 6883148e-04 


Sigma 


xl 


4 


410367359066672e+00 


Sigma 


xdot 


5 


242426264353448e+00 


Sigma 


yi 


6 


995935844827603e+00 


Sigma 


ydot 


7 


681655233650139e+00 


Sigma 


x2 


8 


648715537479415e+00 


Sigma 


par 


8 


795891950965153e+00 


Sigma 


y2 


1 


337220636606977e+01 



Correlation Matrix * 100 



Al 100 
A2 
xl 
xdot 

yi 

ydot 
x2 
par 

y2 

Chi-Sguare 



-61 

100 



1 

-2 
100 



5 -14 
-5 17 
3 
100 -10 
100 



x2 


par 


y2 


2 


-2 


-19 





3 


16 


59 


-13 


14 



8 ii ; 

-18 16 
100 -10 
100 



3255.983023208732 



-10 
1 59 
-3 -1 
-13 -1 
100 3 
100 

DOF = 1246 



Fig. 9. Part of the GaussFit output (slightly edited) ob- 
tained while fitting the double star model in Fig. || to the 
TD for HIP 97237. 



position of the secondary relative to the reference point 
(x2 = Aa*2, yl = A62); and the intensities of the com- 
ponents, Al = A±, A2 = A2. The components are assumed 
to have the same parallax and proper motion. The expres- 
sions within the export () functions are easily recognized 
as the equations of condition, Eq. ([l4]), written in terms of 
the model parameters. The five export () statements are 
divided among two import () loops (which means that the 
data file is forced to be read twice in each iteration): the 
reason is that GaussFit in its standard distribution ver- 
sion cannot handle more than four simultaneous equations 
of condition. 

The model in Fig. || was applied to the TD of 
HIP 97237, using as starting approximation (3600, 3600, 
0, 0, 0, 0.04, 4000, 4300, 0.02) for the variables in the 
parameter list (cf. Sect. |4.3| ). The 'fair' metric with an 
asymptotic relative efficiency of 0.95 was used for robust 



estimation of the parameters (Jefferys et al. 1988a). Part 
of the output file, containing the results of the final (10th) 
iteration, is shown in Fig. ||. It should be noted that the 
estimated standard errors (sigma values) given in the out- 
put file have already been scaled by {x 2 /v) 1 ^ 2 , using the 
chi-square (% 2 ) and degrees of freedom (y) given at the end 
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of the file. Adding the results of the model fitting to the 
reference point data (Sect. 4.3) and using the magnitude 
conversion formula Hp = — 2.5log(A/K) we obtain the 
following estimated parameters of the binary HIP 97237 
(ICRS, epoch J1991.25): 

ax = 296.43961803 deg± 4.41 mas, 

a 2 = 296.43974889 deg ± 8.65 mas , 

8i = 27.12836895 deg ± 7.00 mas, 

S 2 = 27.12856546 deg ± 13.37 mas, 

7T = 80.74 ±8.80 mas 

= -74.94 ± 5.24 mas yr _1 , 

US = -1203.93 ± 7.68 mas yr" 1 , 

H Pl = 12.88 ±0.01 mag, 

Hp 2 = 13.60 ±0.02 mag. 

These data are in reasonable agreement with the values 



derived by Soderhjelm (1999) in an orbital solution com- 
bining the TD with ground-based speckle observations. 



6. Software availability 

Fortran programs interfacing the TD with aperture 
synthesis software and with GaussFit are avail- 
able via the Lund Observatory Internet address 

The 



http: / / www.astro.lu.se / ^lennart /TP / . The program 
td2uv.f extracts the TD for a given HIP identifier and 
converts them into a UV-FITS file that can be used 
e.g. by Difmap. The program td2gf.f similarly extracts 
TD data and generates a data file suitable as input for 
GaussFit. Sample data files, additional information on 
the TD (including descriptions of the known errors), 
and links for retrieving aperture synthesis software and 
GaussFit are also given at this site. 



7. Conclusions 

The Hipparcos Transit Data, published as part of the Hip- 



parcos and Tycho Catalogues (ESA 1997), provide data 



from an intermediate step in the data reduction process 
of the NDAC data analysis consortium. Transit Data are 
given for all known or suspected double or multiple star 
systems in the Hipparcos Catalogue, or about a third of 
the objects in the catalogue. There were several reasons 
to include these data in the published catalogue. A main 
reason was the realization that every stellar system could 
not be examined for every possible type of solution within 
the time available for completing the catalogue. Ideally, 
such examination should also take into account ground- 
based observations. The Transit Data allow the user to 
re-examine such solutions, should the need arise. 

Basically the Transit Data contain the Fourier coeffi- 
cients which describe the modulation of the detector signal 
caused by the object's motion across the modulating grid. 



As such they retain all the photometric and astrometric in- 
formation on the object gathered during its transit across 
the grid. The Transit Data have been carefully calibrated 
and referred to the Hipparcos photometric and astromet- 
ric systems, so that any data derived from them should 
be directly comparable with other results in the published 
catalogue. 

By giving this review of how the Transit Data were 
recorded, what they physically represent and examples of 
their practical uses, we hope to encourage readers to uti- 
lize the data in their own exploitations of the Hipparcos 
results. To aid this process, we have made programs avail- 
able which provide interfaces with publicly available soft- 
ware packages, in particular Difmap and GaussFit. Many 
other applications could be thought of - the combination 
of Transit Data with ground-based speckle observations of 
double stars to improve orbits, parallaxes and mass ratios 



by Soderhjelm (1999) is an example. What has been cov- 
ered and demonstrated here might just be a stepping stone 
to new and creative exploitations of the Transit Data. 
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